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1 Introduction 



Of late the rare decays of the B mesons have been recognized as important 
tools to study the basic structure and validity of the Standard Model (SM) 
and its extensions. In particular, the radiative decays, owing to their relative 
cleanliness as far as experimental signatures are concerned, have attracted a 
great deal of attention. For a general overview of the kind of issues considered 
relating to radiative decay modes, see [H El EI and references therein. The 
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decays B — > X^'j and B K*'j have been observed |2] and both these 
are extensively used and rehed upon for constrining the parameters of any 
new theory or extension of the SM The decay B —>■ K*'-/'-/ is another 
potential testing ground for the effective quark level Hamiltonian b — 577 
first studied by Lin, Liu and Yao |1] and pursued further in references |S],|n] 
and [7] . This amplitude has been the focus of considerable research recently, 
not only for the useful indications it will give to the underlying theories of 
flavour changing neutral currents, or the possible contributions from loops 
with super symmetric partner particles, but for the impending experimental 
studies of the B-factories in the near future. 

As has been previously noted, the b 577 amplitude naturally splits into 
two categories: an irreducible contribution which is well known and usually 
estimated through basic triangle graphs, and a reducible one, where the sec- 
ond photon is attached to the external quark lines of the b —>■ sj amplitude. 
At the quark level the reducible contribution presents no real problems, how- 
ever, when we consider an exclusive channel, such as i? ^ M77 for a specific 
meson M, it becomes more appropriate to consider the second photon as 
arising from the external hadron legs of the amplitude B —>■ M7. In contrast 
to the earlier cases of B ^ i^(7r)77 |7] where the amplitude for a single 
real photon vanishes identically, resulting in the irreducible diagram to be 
the sole contributor, the amplitude B K*'j is non vanishing. However for 
the neutral decay mode the second photon cannot arise from the resulting 
K*, and thus in this case also, it is the irreducible amplitude that stands 
out, though for a completely different reason. Of course we must also con- 
sider the usual long distance contributions, such as the process B —>■ K*r]c 
followed by the decay rjc 77. Note that for completeness we will also in- 
clude the rj contribution, even though the 1] coupling to cc will be small. The 
rate for B — »• Kt]' is anomalously high and many possible mechanisms have 
been proposed that aim at taking this anomalous production into account |8j. 
However, in the present case there is not enough data corresponding to the 
B K*!]' channel, and at present only an upper limit on this branching 
ratio is available. We therefore tend to remain conservative in the present 
study regarding this issue and assume that the rj' contribution can be ob- 
tained similar to the rjc contribution. The situation is expected to improve 
with the availability of more and precise data in this direction. We therefore 
include an rj' contribution along the lines of the r]c contribution. 

In this paper we will estimate the branching ratio for the process B^ — >■ 
i^'*°77 by considering the effects of the irreducible triangle diagram contribu- 
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tions in the next section, followed by the resonance contributions in section 
3. Note that in the case of -B — 7^*77 there will only be three sizeable res- 
onance contributions; r/c, f] and r]' . Furthermore, each of these contributions 
will only contribute a narrow peak in the 77 invariant mass spectrum, which 
is easily separated experimentally. As such the interference terms for each of 
these pairs of terms will not be considered here. Finally in section 4 we will 
present our results and analysis. 



2 The Irreducible Contributions 



The irreducible triangle contributions to the process in which we are in- 
terested [B — >• iir*77) originate from the quark level process h S'j'j. The 
effective Hamiltonian for this process is J3] 
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The invariant amplitude corresponding to this effective Hamiltonian is 
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where 



and 



T^u = k2p,kiy — {ki.k2)gi_iu, 
Au = 3(^3 - C,) + (C4 - a); Ad = \Au, 
A, = 3(Ci + C3 - C5) + (C2 + - Ce), 
As = A, = i[3(C3 + - C5) + (C3 + C4 - a)], 
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i? = C = --(3C6 + C5). 



(4) 
(5) 



(6) 



In the above expressions we introduced the functions 

jK) = JnK), irK) = 4/nK)-/ooK), L^) = /qoK), 
where 



1 rl-a; 



, • (7) 

70 m? — 2{ki.k2)xy — ie 

Note that to get the M.{B K*'-/'-/) invariant amphtude from the quark 
level amplitude we replace the (s|r|6) by (i^'*|r|i?) for any Dirac bilinear F. 

With q = Pb — Pk* = ^1 + ^2 and following Cheng et al. we parame- 
terize the hadronic matrix elements as 



{K*{pK^)\-s^,h\B{pB)) 
{K*{pK^)\s^,l,h\B{pB)) 



e,ua(3e*''{pK')p"BP^K', (§) 



mB + rriK 
{mB + mK*)e*^'{pK*)Ai{q^) 

— {Pb +PK*)pA2{q'^) 

niB + niK* 

mK 



(9) 



-2^{,\pB)q,{A,{q^)-A,{q')} 

For the functional dependence of various form factors appearing above, we 
follow PO]. Using these definitions, we determine the irreducible matrix 
element for the process B 7^*77 as 
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where 



pup 
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The functions Ki above are defined as 
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3 Resonance contributions 

For this process there will be three significant resonance contributions, 
that from the r]^--, r]- and T^'-resonances. The rj^ contribution to the decay 
process comes via the t-channel decay B — > K*ric, with the rjc then decaying 
into two photons. 

The T-matrix element for this process can be written as 



(ir*77|r|5) 



{K*r^,\T\B){^^n,) 



g2 - + im^Xtotai 
The amplitude (77|T|?7c) is parameterized as [U] 

(77|T|r^.) = 2iB,X''"Pel{k,)el{k2)k,^k2p. 



(12) 
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Note that we can determine -B^^ from the known decay rate: 



(27r)32A;? (27r)32A;^ 

(14) 



where we have 



and so 



Y: \{i^n.)f = 2\B,fq\ (15) 



spins 



- 77) = (16) 
lore 

The -B ^ K*r]c amphtude has been determined in Cheng et al. [9^ as 



M{B K*vc) = ^y^V:, ( C, + -C, ) '^^> (17) 



where 



^(Boi.-,..) ^ 2/,,m;,*<^*(my(e;,,.pB). (18) 

Note that /^^ is defined as (0|c7^c|r7c) = ifr^.PfiiVc)- 

In the above way of parametrizing the B K^rjc, there is a lot of model 
dependence that goe in. Since, the branching fractions for this sub-process 
is known, we can in principle avoid such a model dependnce by writing the 
amplitude as 

M{B^K*r^,) = aff^{e*^,.pB) (19) 

and determine the effective constant from the corresponding decay rate. We 
folow this procedure and therefore try to avoid any model dependence as far 
as possible. 

Therefore the total contribution due to the rjc resonances is thus, 

= ^B,^^Tff . jrr^^^ r''- e^-'-'el^elMahp. (20) 

y — ^r)c "''■r^c-'- total 

Analogous to the rjc resonance, the 1]- and T^'-resonance contributions, Airj 
and Air)', have exactly the same form as equation pn|l with the parameters 
Brj^, m.^^ and T'^^^^i being replaced by their rj- and r^'-counterparts respectively. 
However in this case we cannot define the relative sign of the amplitudes and 
any of the other components of the amplitude. 

The process can also receive additional contribution from the B* and K2 
channels where the B meson decays into a photon and an on-shell K2 or 



6 



slightly off-shell B* and then these giving rise to K* and the second photon. 
However, there is no data available at present for either of these and if the 
widths for the individual channels contributing to the process are significant, 
the contribution can be sizeable. However, we expect that these contributions 
can be eliminated by suitable cuts in the B*- (or K2-) photon plane and thus 
we do not consider them here at all. 



4 Results 



The squared amplitude for the process B — K*jj is then; 
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where the interference terms have not been included here. The components 
to the squared amplitude were calculated to be; 
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AVe have similar expressions for the rj and rj' terms, replacing the parameters 
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by their r] and 7]' counterparts. 
The irreducible squared matrix element is then; 
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The total decay rate is then given by 



dV 




1 1/2 



d{cos9)dy^. 
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(26) 



where -y/s^ is the CM. energy of the two photons while 6 is the angle which 
the decaying B-meson makes with the two photons in the 77 CM. frame. 

Our results arc presented in figures 1 and 2 (both plotted with a loga- 
rithmic scale on the y-axis), where figure 1 shows the differential branching 
ratio given as a function of the invariant mass of the two photons, for the 
case of the neutral B meson decay without the interference terms between 
the resonances and the irreducible background included. The inclusion of the 
interference terms will in principle give rise to an interference pattern near 
the base of the resonance peaks. Since the peaks are narrow and moreover 
the interference contributions being small (except the r^c-Irreducible term), 
we do not include them in the plots. It is worth mentioning that imposition 
of suitable cuts in the spectrum to eliminate the resonance contributions will 
eliminate any such interference patterns also. The numerical estimate for 
the branching ratio arising due to all possible contributions is summarized 
in Table 1. Quite evidently, the largest contribution comes from the 77 reso- 
nance mode. It should be stressed again that using appropriate cuts in the 
spectrum, the resonances can be completely eliminated and what is left is 
the background irreducible contribution. In estimating the numerical values 
for the interference terms, we have assumed that the relative signs between 
the terms are such that the 77-Irreducible and Ty'-Irreducible interfering con- 
tributions add on to the other pieces. However, because of the smallness of 
these values, it really makes no significant difference. 

In figure 2, we compare the irreducible contribution with the total con- 
tribution to the branching fraction. Clearly, the resonances dominate the 
results. At this point, it may be worth mentioning that a quick look at the in- 
dividual values tabulated in Tablel reveal the following. Since the resonances 
are narrow, one may try to estimate the contributions directly by multiplying 
the individual branching ratios ie we expect in the narrow width approxima- 
tion that BR{B > /s:*77) - BR{B - - > K*X)BR{X > 77) 

where X denotes any of the resonances. The numbers quoted clearly show 
that they are in accord with the expectations. However, if we had used the 
form as in Eq(17), we might have over- or under-estimated the resonance 
contributions (except probably for the r] mode) because it is not very clear 
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Contribution 


Branching ratio xlO ^ 


Resonance 




Vc 


4.7 


V 


56.9 




3.7 


Irreducible 


2.3 X 10-2 


Interference 




?7c-Irreducible 


2.6 


?7-Irreducible 


3 X 10-3 


r^'-Irreducible 


4.5 X 10-3 


BR 


57.5 



Table 1: Contribution to the S° — > ir*°77 branching ratio. For the interfer- 
ence terms we quote the absolute values. 

if such a simple parametrization is the correct one. We however avoid any 
such possible conflict by drawing heavily on the experimental values (upper 
limit for t]') of the various sub-process branching fractions. 

The central values of the parameters used in our calculation are shown 
in the Appendix. We make an attempt to estimate the errors creeping into 
the numerical calculations due to errors in various input parameters. The 
theoretical uncertainties arising out of uncertainties in the parameters are 
overwhelmingly in the input values of the form factors and the meson decay 
constants . The Wilson coefficients have been taken to be their NNL values 
and there are no significant theoretical uncertainties in them. In evaluating 
the uncertainties of our results, it is appropriate to evaluate them separately 
for the background irreducible contribution and the resonance contributions, 
since the latter can easily be experimentally separated from the former by 
suitable cuts in the spectrum. For the model dependent parametrization of 
Eq(17), the theoretical uncertainty in the resonance contribution due to rjc 
arises mostly because of uncertainties in the CKM parameters, /^^, B^^ and 
Fo(m^2); the Wilson coefficients values used are the NNL level values and 
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no comparable uncertainties exist therein. The CKM parameters relevant 
to us have an uncertainty of about 10% jTTj. The form factors used are 
the same as in |9| where the actual dependence of the form factors as a 
function of momentum transfer squared are given and hence no errors arise 
due to parametrization of form factors as a function of . Although this 
reference does not quote any estimate of the uncertainties in the numbers, a 
typical uncertainty in this type of calculation based on quark model is given 
in [12] and is typically of the order of 15% , arising to a great extent due 
to uncertainty in the strange quark mass. The rate of the rjc decay into two 
photons is uncertain by about 40% [TTJ . A typical estimate of the uncertainty 
in the value of /^^( arising mostly again out of uncertainty in the current mass 
of the s-quark ) has been estimated at about 15% in [T3]. Combining all this 
, we would expect that an estimate of the contribution of the rjc resonance 
based on such a parametrization to be uncertain by about 50%. A similar 
estimate for the other two resonaces give a somewhat lower value mostly 
because their decay rates into two photons is better known, to an accuracy 
of about 10 % for rj' and about 5 % for the r]. The uncertainties in the decay 
constansts of the eta' and the t] have been estimated to be about 10% [13] 
and we estimate the overall uncertainty in our calculation for the rj' and t] 
to be about 40% and 30% respectively. However, since we have relied on 
experimental values of rates and branching fractions, the above mentioned 
uncertainties are significantly reduced and the only source of uncertainty in 
our estimation is the uncertainty present in sub-process rates. 

Turning to the irreducible contribution, the uncertainty arises mostly 
because of the CKM factors and the form factors. These combine to give an 
overall uncertainty of about 20% for the irreducible part of the amplitude. 
As stressed before, once the suitable cuts are imposed in the two photon 
spectrum, it is possible to extract the irreducible contribution and here the 
errors are relatively smaller and are expected to even go down further with 
more accurate determination of the CKM parameters and the form factors. 

At the levels reached by the current B-factories, the branching ratios 
obtained are too low to be observed. One certainly hopes that in the near 
future experiments, with better luminosities possible, the numbers obtained 
will be very useful for confronting theoretical models with experimental data. 
As discussed in the text, this decay with two photons depends on the parts 
of the effective Hamiltonian, which the decays with a single photon are not 
sensitive to and thus provides a more complete test of the underlying theory. 
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5 Appendix 

We list the central values of the various parameters entering our numerical 
estimates: 

GF^l.lGxlQ-^Gev-^ 0(7715) = 1/130 mb = 4:.8Gev 
rrit = 175Gev rric = 1.5Gev rriu = = 777^ 
mK,+ = 0.89Gev mK*o = 0.896^67; 

F(f ^* = 0.3 Mpoie = 6.65^67; 
FtotiVc) = 1-3 X 10-2 rur,, = 3.06-67; 
Br,^ = 2.74 X 10-^^67;-^ fr,, = 0.35^67; 
Ttotiv') = 0.203 X 10-^^67; m^/ = 0.96^67; 
Br,' = 14.0 X 10"-'^G'e7;-^ = -6.3Me7; 

B{r]' 27) = 2.11% B{r]c ^ 27) = 3 x 10"^ 
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Figure 1: The differential branching ratio oi B ^ 7^*77 plotted as a function 
of the CM energy of the diphoton rest frame without interference terms taken 
into consideration. Here we have plotted with a log scale on the y-axis. 
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Figure 2: The irreducible and the total contribution to the differential 
branching ratio of i? ^ 1^*77 plotted as a function of the CM energy of 
the diphoton rest frame. 
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